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ABSTRACT 

We examine the dynamical evolution and statistical properties of the supernova ejecta of massive 
primordial stars in a cosmological framework to determine whether this first population of stars could 
have enriched the universe to the levels and dispersions seen by the most recent observations of the 
Lyman-a forest. We evolve a ACDM model in a 1 Mpc'^ volume to a redshift of z = 15 and add 
"bubbles" of metal corresponding to the supernova ejecta of the first generation of massive stars 
in all dark matter halos with masses greater than 5 x IO^Mq. These initial conditions are then 
evolved to z = 3 and the distribution and levels of metals are compared to observations. In the 
absence of further star formation the primordial metal is initially contained in halos and filaments. 
Photoevaporation of metal-enriched gas due to the metagalactic ultraviolet background radiation at 
the epoch of reionization (z ~ 6) causes a sharp increase of the metal volume filling factor. At z = 3, 
~ 2.5% of the simulation volume (« 20% of the total gas mass) is filled with gas enriched above a 
metallicity of IQ^^Zq, and less than 0.6% of the volume is enriched above a metallicity of 1Q~^Zq. 
This suggests that, even with the most optimistic prescription for placement of primordial supernova 
and the amount of metals produced by each supernova, this population of stars cannot entirely be 
responsible for the enrichment of the Lyman-a forest to the levels and dispersions seen by current 
observations unless we have severely underestimated the duration of the Pop HI epoch. However, 
comparison to observations using carbon as a tracer of metals shows that Pop III supernovae can be 
significant contributors to the very low overdensity Lyman-a forest. 
Subject headings: cosmology: theory ~ intergalactic medium - galaxies: formation 



1. INTRODUCTION 

Recent observations by iSchave et all l)2003(l have 
shown that the Lyman-a forest is polluted with metals 
at very low densities. The distribution of metal is very 
strongly dependent on overdensity, with median metal- 
licity values ranging from [C / H] — —4.0 at log 5 — —0.5 
(where 5 is defined as [5 = pjp) to [C/i?] — —2.5 at 
log 8 = 2.0 using their fiducial UV background model. 
Their observations show little evidence for metallicity 
evolution of the Lyman-a forest over the redshift range 
z = 1.5-4.5. 

The lack of observed evolution in metallicity is sugges- 
tive of a very early epoch of stellar evolution. Recent ob- 
servations by the Wilkinson Microwave Anisotropy Probe 
suggest an epoch of star formation in the redshift range 
of z = 11 — 30 (,Kogut et al. .2003.), which is co nsis- 
tent with the s imula tion results of lAbel et alJ l)2002H and 
iBromm et al.l ()2002|) . which suggest that the first gen- 
eration of stars (known as Population HI, or Pop III) 
formed in the redshift range z = 20 — 30. The Abel 
et al. results, which are the highest-resolution simula- 
tions of formation of the first generation of primordial 
stars to date, also suggest that Pop HI stars are very 
massive - on the order of ^ 200 A/©. Stars that are 
in this mass range will die in extremely energetic pair- 
instability supernovae and can eject up to 5 7 M© of ^^Ni. 
l|Heger et al.ll2?)03l iHeger fc Wooslevll20?)^ . The forma- 
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tion site of Pop III stars is in halos with total m asses of 
~ IO^Mq. l|Abe l et al.'2002rYos hida et al .'2003). These 
halos have escape velocities which are on the order of 1 
km/s. Due to the sh allowness of the potential wells that 
Pop III stars form in. iFerraral l)1998t) suggests that ejecta 
from a massive Pop HI supernova can propagate to very 
large distances (far greater than the halo virial radius), 
a result wh ich is supported in simulations performed by 
IBromm et a l. (2003). 

In this paper we describe the results of cosmological 
hydrodynamic simulations which address whether or not 
a population of massive primordial stars can be respon- 
sible for metal enrichment of the Lyman-a forest to the 
level and dispersion seen today. We examine the most op- 
timistic possible scenario for Pop III star formation and 
enrichment in order to establish an upper limit on metal 
enrichment of the Lyman-a forest due to Population III 
stars. 

2. simulations 

The simulations were set up using the concordance cos- 
mological model {VL,n ~ 0.3,S1a = 0.7, fif, = 0.04 and a 
Hubble parameter oih = 0.7 where h is in units of 100 km 
s~^ Mpc^"'^. Initial perturbations in the dark matter and 
baryon density were created with an lEisenstein fc Hul 
(,1999) power spectrum with (Tg = 0.9 and n = 1 us- 
ing t he Zel'dovich approximation ijBertschinger fc Gelbl 
|1991|) in a box which is 0.7 comoving Mpc on a side. 
In our simulations we use a computational box with 256'^ 
grid cells with a comoving spatial resolution of 2.7h~^ 
kpc and a dark matter mass resolution of 1477 Mq. 

The simulation was initialized at z=99 and was allowed 
to evolve to z=15 us ing the Euler ian adaptive mesh re- 
finement code Enzo l)BrvarJll996j) . The simulation was 
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stopped at z — 15 and all halos with dark matter mass 
M_DM > 5 X 10^ M(7) were found using t he HOP halo- 
finding algorithm ijEisenstein fc Hutlll998D . This yielded 
184 halos, each of which we assume produces one pair in- 
stability supernova. Note that we are ignoring negative 
feedback, which rais es the minimum halo m ass which can 
form Pop III stars ijMachacek et al.ll2001j) . This is con- 
sistent with our desire to simulate a best-case scenario. 
We discuss this and other assumptions in Section 0] 

At this point, two separate cases are considered. In 
the first case (referred to as Case 'A'), spheres of uni- 
form metal density 1 kpc (proper) in radius with 127 
Mq of metal are placed in the simulation, centered on 
the cell of highest dark matter density of each halo. This 
corresponds to a mass averaged metallicity in the volume 
oi<Z >= Mz/Mb = 4.02 X 10"^ Zq (where Mz and 
Mb are total metal and baryon masses in the simulation 
volume), which remains constant throughout the simula- 
tion. No other modifications to the data were made - in 
particular, the baryon density, temperature and peculiar 
velocities and dark matter density were unmodified. 

In the second case (Case 'B'), the spheres of uniform 
metal density are placed down in the same manner. In 
addition, the baryon gas in the corresponding volume is 
smoothed to the cosmic mean (< p >= f^bPc), and the 
temperature of the baryon gas is set to 10"* K. This cor- 
responds to the net smoothing and heating of baryons in 
primordial halos due to pair-instability supernovae. The 
mass averaged metallicity in this case is slightly greater, 
< Z >= 4.11 X 10~^ Zq. This is due to a small net loss 
in baryon mass when the densities are smoothed. As in 
Case A, peculiar velocities and dark matter density are 
unmodified. 

The simulations are then evolved to z=3, following 
the evolution of the gas using non-equilibrium chemistry 
and cooling from hydrogen and helium and their ions 
ifAbel et al.lll997HAnninos et alJll997|) . Metal line cool- 
ing is ignored. At z = 7 we initializ e a uniform meta- 
galac tic ultraviolet background field ijHaardt k. Madaul 
11996^ with a slope of = —1-8 that is ram ped to a 
pre-c omputed value by a redshift of z = 6. Zha ng et al.l 
()1997l|) determined that such a prescription can reproduce 
the Gunn-Peterson effect in the hydrogen Ly-a lines by 
a redshift of 5.5. Self-shielding of the gas is not explicitly 
included in the assumed UV background. 

It is important to note that the box size, at 0.7 h~^ 
Mpc, is somewhat small. Statistical results are reliable 
at z = 15. However, by z = 3 (when the simulations 
are terminated) the box is too small to be statistically 
reliable. We also performed a simulation in eight times 
the volume at the same mass and spatial resolution (i.e. 
512^ cells/particles), which gave results indistinguishable 
from what follows. Nonetheless, all results at late times 
should be considered qualitative. 

3. RESULTS 

The general evolution of our simulation after injection 
of metal bubbles at z = 15 is as follows: Beginning at 
z = 15, the bubbles of metal track the flow of gas onto 
and along dark matter filaments. The competing effects 
of advection along the filaments and the collapse of fil- 
aments during this period essentially cancel out, with 
little net effect on the fraction of the volume occupied 
by metal enriched gas (also referred to as the volume 



filling factor, or VFF). Regions of relatively high metal- 
hcity [Z > IO^'^Zq), corresponding to the densest re- 
gions of filaments, decrease their volume filling factor 
significantly from z=15 to z=6. Case A shows more of a 
decrease in VFF at relatively high metallicities (from a 
VFF of 10-2-5 to 10-3-2) than Case B does (which has a 
minimum VFF of 10-^-^) due to the higher initial den- 
sities and lower initial temperatures of polluted regions 
in Case A. Figure ^ shows snapshots of the metal dis- 
tribution in the simulation volume taken at 4 different 
redshifts, and Figure [21 shows the volume filling factor 
for three metallicity thresholds as a function of redshift. 

A uniform metagalactic ultraviolet background is 
switched on at z = 7. Photoheating raises the mean tem- 
perature of the baryon gas. In the Lyman-a forest over- 
density regime (1 < S < 10), which roughly corresponds 
to the filamentary structure observed in Figure ^ the 
temperature is raised to ~ 10,000 K. The local thermal 
speed of the baryon gas then exceeds the escape veloc- 
ity of the filaments, resulting in significant expansion of 
the volume occupied by the gas in those filaments. This 
includes the gas polluted by metals. This effect can be 
clearly seen in the third and fourth panels of Figure ^ 
Figure|21shows the sharp increase in metal VFF in a more 
quantitative way. Gas which has been polluted above a 
metallicity of 10-^^0 (corresponding to essentially the 
entire volume of gas polluted by metals) increases in VFF 
to 0.048 for both Case A and 0.046 for Case B. This 
corresponds to 28% (26%) of the baryon mass being en- 
riched for Case A (Case B). The two cases are essentially 
indistinguishable with regards to the total mass and vol- 
ume of gas polluted by metals. Examination of gas with 
higher metallicity {Z > 10~^Zq) shows some difference 
between the cases, with a maximum VFF of 5.2 x 10^"^ 
for Case A and 5.8 x 10--^ for Case B (corresponding to 
1.5% and 1.9% of the total gas mass, respectively). This 
difference is due to the initial smoothing and heating of 
baryons in Case B. 

In Figure |2| we estimate the amount of carbon con- 
tained within the primordial metallicity field at z = 3. 
We assume that the carbon abundance in the metal den- 
sity is equal to Xq = 0.027, which is taken from the 
supernovae metal yield of a massive pri mordial star of ap- 
proxim ately 260 Mq as computed bv iHeger fc WooslevI 
l)2002j) . We then compute the mean and median car- 
bon metallity, [C/H] = log {nc/riH) — [C/H]q (where 
[C/H]q — —3.45) in bins of constant logarithmic over- 
density between log S = —0.5 — 2.0 and plot this in 
Panel A of Figure |3| Altering the effects of Xc results 
in this figure being scaled along the y-axis by a factor 
of log(Xr7/0.02 7). The solid green line is the fit to the 
observations of iSchave et alJ l)2003(l using their fiducial 
model. The dashed-green line is the fit of the lower bound 
lognormal scatter of their data. The results of our sim- 
ulation yield that the Population III carbon content in 
the IGM at z = 3 is below the observed limits across the 
entire overdensity range. 

In Panel B of Figure|2|we plot the probability distribu- 
tion function (PDF) for both simulation cases in the over- 
density range log S = —0.5 to -1-2.0. The vertical lines 
correspond to the mean (solid) and median (dashed) val- 
ues of the [C/H] values. The distributions corresponding 
to the two cases are statistically indistinguishable within 
one standard deviation. The small variation between the 



Fig. 1. — Projected log metal density (Case A). The area viewed is a projection of the entire simulation volume. The four panels 
correspond to (from left to right) z=15, 7, 5 and 3. The Pop III supernova remnants are placed in the volume at z=15 and advect along 
the filaments. Photoevaporation of gas in the filaments, driven by the metagalactic UV background, causes the volume filling factor of the 
metal-enriched gas to increase substantially by z=3. 
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Fig. 2. — Volume filling factor. The lines describe the fraction of the simulation volume filled to a metallicity of at least 10~® Zq (black), 
10~^ Zq (blue) and 10~^ Zq (red). Panel (a) corresponds to the simulation where spheres of uniform metal density are added and no 
other changes are made. Panel (b) corresponds to the simulation where, in addition to uniform spheres of metal density, the baryon density 
in the corresponding volume is smoothed to the mean density of the simulation and the gas temperature of the sphere is raised to 10* K. 



two cases in the range [C/H]=-4 to -2 is due to the differ- 
ence in their mean initiai metaliicity per bubbie (caused 
by the difference in treatment of of the baryon density 
ficid in the region initiaiUy poUuted by nietais in the two 
cases) . 

In order to determine the observabiiity of the primor- 
dial metallicity field, we postprocess our data to com- 
pute the fraction of CIV to neutral hydrogen (HI) for 
each cell in our computational volume. Obtaining the 
quantity log F — log (jiqiv /'tihi) within constant over- 
density bins, in the range log 5 = —0.5 — 2.0, allows the 
determination of the lognormal average < log F > at 
each redshift. At z=3 our analysis yields < log F >= 
-3.78 ± 0.91(-3.78 ± 0.92) for Case A (B). We can 
then approximate the CIV optical depth due to PopIII 
stars as r^fy^^ « ( '^7^x17 ) ' 10"^'°^ ^^thi, where 
fciv.Hi ^civ.Hi are the oscillator strengths and 
rest-frame absorption wavelengths for the two species 
(CIV & HI). Using our mean values for < log F > at 
z=3 we obtain an estimated optical depth, due to CIV, 
of r^rr/" = 0.971!^ X 10-4 THi {0.97l:ll_x 10"^^/). 
iSchave et alJ l|2003l') computed tciv = thi in the 

Ly-a forest overdensity range. If the above value corre- 
sponds to the total metallicity at z=3 then the contri- 
bution of the primordial component to the total optical 
depth of CIV is about r^"/^" = 0.01 tciv- This resuh is 



some what sensitive to the shape of the UV background 
- see ISchave et abl ll200^ for more details. Statistical 
correlations between PopIII CIV and HI absorbers and 
more detailed examination of spectra due to ejecta from 
Population HI stars will be discussed in a forthcoming 
paper. 

4. DISCUSSION AND CONCLUSIONS 

In this paper we use cosmological hydrodynamic simu- 
lations to examine the evolution of metals ejected by an 
early population of massive primordial stars. We show 
that, in the absence of further star formation, photoe- 
vaporation of baryons bound to dark matter filaments 
during rcionization is the most important mechanism in 
determining the volume filling fraction by z = 3. Our 
two study cases, although different in their initial setup, 
give the same results for the global distribution of the 
primordial metal field by 2; = 3, suggesting that our re- 
sult is insensitive to small-scale dynamics. 

Comparison of our r esults to obser vations of carbon in 
the Lyman-a forest bv ISchave et al.l lf2003) show that at 
z = 3 the median value of the Pop III carbon metal- 
licity for both cases considered fall within the low end 
of the scatter range of the observed data for log 6 < Q. 
For log (5 > the Population III carbon metallicity is 
below the observed values, with the Schaye result show- 
ing a much stronger increase in metallicity with overden- 
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Fig. 3. — In panel (a) we plot the volume averaged [C/H] in solar units within constant logarithmic overdensity bins in the range of 
1 < (5 < 10^ at z = 3. The profiles for the two cases discussed in the text show no statistical difference. In panel (b) we plot the probability 
distibution function (PDF) of [C/H] within the overdensity range 1 < <5 < lO^. The y-axis measures the fraction of the metal polluted 
volume with [C/H] values between [C/H] and [C/H] -|- A [C/H]. The result is computed with Xc = 0.027. The two distributions yield 
< [C/H] >= -4.68 ± 0.80 for Case A and < [C/H] >= -4.68 ± 0.81 for Case B. Median values for the two PDFs are -4.53 and -4.56 
respectively. The two cases are indistinguishable within one standard deviation. 



sity, resulting in the median value of [C/ H]^°p^^^ becom- 
ing an increasingly smaller fraction of the total observed 
[C/H]. 

Our results depend strongly on two factors, namely, 
the total number of Population III stars formed in our 
volume and the metal yield per star. In these simula- 
tions we make the assumption that all halos with mass 
Mum > 5 X IO^Mq form a massive primordial star by 
z = 15, which was gu ided by the simulat i ons pe rformed 
bv lAbel et all l|2n02D and lYoshida et alJ l|2003D . which 
show that this is the characteristic dark matter mass of 
a halo w hich forms a star in th e early universe. Simula- 
tions by 'Machacek et al. (2001) and semianalytical cal- 
culations by Wise fc Abel (2003) show that a soft UV 
background produced by the first Pop III stars effectively 
dissociates H2, which is the primary cooling mechanism 
in primordial star formation. This so-called negative 
feedback effect raises the minimum halo mass that can 
form a primordial star within it and therefore reduces 
the number of ha los which will form P opulation III stars 
at a given epoch. iWise fc Abell ()2003D find that negative 
feedback reduces the number of star forming halos by a 
factor of 5-10 relative to what we used. On the other 
hand, suppression of Pop III star formation by negative 
feedback would be compensated by an extended epoch 
of Pop III star formation. At present, we do not know 
when Pop III star formation ceases. We view our choice 
of Mrmn = 5 X lO^M© at z = 15 as a hedge between 
competing effects. 

The decision to place spheres of metal in the simulation 
volume at 2; = 15 was guided prirn arily by the WMAP 
polarization results. (jKogut et al.l 12003). This choice 
may have resulted in an underestimation of the number of 
Population III stars (and therefore metal pollution due to 
Pop III supernovae) because there are dark matter halos 
which form after z = 15 b ut may stil l be unpolluted by 
metals. However, results bv lBromm e t al. ( 2001) suggest 
the existence of a "critical metallicity" of 5 x 10~^Zq 
above which a solar IMF dominates, and it has been 
argued that this metallicity is reached by z ~ 15 — 20. 



(jMackev et al.ll2003l:ISchneider et a'l]l2002(l The choice of 
z = 15 for our epoch of instantaneous metal enrichment 
seems to be a reasonable compromise. 

The physical properties of the metal "bubbles" can 
have a possible effect on our results. The choice of a 1 kpc 
(proper) radius for the metal bubbles is somewhat arbi- 
trary. Several calculations have been performed that sug- 
gest that ejecta from the most massive pair -instability su- 
pernovae can propagate to large distances ijBromm et alJ 
l2003tlMadau et al.l200lD . but the n iaximum propagation 
distance is unclear. Additionally, iBromm 
suggest that the ejecta from pair-instability supernovae 
still has substantial peculiar velocities (~ 50 km/s) at 
500 pc. The metal spheres in this calculation have no 
initial outward peculiar velocity, which may result in a 
smaller volume filling factor than if this were taken into 
account. 

The second factor that strongly affects our result is 
the choice of the a mount of m e tals c reate d per Popula- 
tion III supernova. lAbel et all (|2002t) and IBromm et alJ 
l|2002i) both suggest that the first population of stars will 
be very massive. The mass function of the first genera- 
tion of stars is unclear, due to lack of resolution and ap- 
propriate physics. The main-sequence mass of the star 
strongly affects its ultimate fate: Stars with the range 
of ^ 140 — 26OM0 detonate in pair instability super- 
novae, which are much more energetic (up to ~ 10^'^ ergs 
compared to 10^^ ergs for a standard Type I or Type 
II supernova) and produce more metal (up to 57 Mq of 
s^Ni and almost 130 Mq of total metals for a 260 Mq 
primordial star). However, stars between ~ 50 — 140Mq 
and above ~ 26OM0 form black holes without first eject- 
ing significant quantities of nucleosynthesized material. 
(Heger ct al. 2003). The amount of metals placed into 
the simulation volume is scalable - if the mean amount 
of metals ejected by Population HI stars were lower (due 
to some substantial fraction collapsing directly into black 
holes, for instance), all of the results shown in Figure [21 
and Panel A of Figure O scale linearly with the mean 
amount of metal produced per star. 
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Our results for [C/H] vs. overdensity (using carbon 
as a proxy f or me tallicity) agree with the results of 
iSchave et alJ l)2003|) to within one standard deviation for 
the lowest observed densities (log (5 < 0). These are the 
densities that are the most likely to remain unpolluted 
by later generations of stars, which form in deeper poten- 
tial wells. Further study of the lower density regions of 
the Lyman-a forest could yield more constraints on the 
mass and total number density of massive Population III 
stars. 

An additional factor to consider is that the nucleosyn- 
thetic yields of very massive primordial sta rs are much 
different than that of metal-enriched stars. l)Heger et al.l 
|2003). Due to this, it may be possible to disentangle 
the effects of massive primordial star s and th eir metal- 
polluted descendants, as discussed bv lOh et"al . (2001). 

Due to our choices of low minimum halo mass and high 
metal yield per supernova, our result is a strong upper 
limit on the pollution of the Lyman-a forest due to Pop- 
ulation III stars, unless we have severely underestimated 
the duration of the Population III epoch. 

The simulation volume is relatively small. Though a 
reasonable statistical representation of the universe at 
z = 15, the results obtained at later times {z ~ 3) 
should be considered qualitative due to the small box 
size. A much larger simulation volume is required for ad- 



equate statistics at z 3. However, simulating a much 
larger volume which would still have reasonable spatial 
and dark matter mass resolution on a single grid is com- 
putationally prohibitive at the present time. 

All of the simulation results described in this paper are 
performed without further star formation or feedback. 
Having a single episode of star formation at z = 15 means 
that metal evolution after that time is passive, whereas 
in reality there would be continuous star formation and 
feedback. A logical extension of this work is the inclusion 
of later epochs of star formation and their resulting feed- 
back of metals and energy into the IGM. These results 
will be presented in a forthcoming paper. 
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